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Haberlea rhodopensis is a resurrection species with extreme resistance to drought stress 
and desiccation but also with ability to withstand low temperatures and freezing stress. In 
order to identify biochemical strategies which contribute to Haberlea's remarkable stress 
tolerance, the metabolic reconfiguration of H. rhodopensis during low temperature (4°C) 
and subsequent return to optimal temperatures (21 °C) was investigated and compared 
with that of the stress tolerant Thellungiella halophyla and the stress sensitive Arabidopsis 
thaliana. Metabolic analysis by GC-MS revealed intrinsic differences in the metabolite 
levels of the three species even at 21 °C. H. rhodopensis had significantly more raffinose, 
melibiose, trehalose, rhamnose, myoinositol, sorbitol, galactinol, erythronate, threonate, 
2-oxoglutarate, citrate, and glycerol than the other two species. A. thaliana had the 
highest levels of putrescine and fumarate, while T. halophila had much higher levels 
of several amino acids, including alanine, asparagine, beta-alanine, histidine, isoleucine, 
phenylalanine, serine, threonine, and valine. In addition, the three species responded 
differently to the low temperature treatment and the subsequent recovery, especially 
with regard to the sugar metabolism. Chilling induced accumulation of maltose in H. 
rhodopensis and raffinose in A. thaliana but the raffinose levels in low temperature 
exposed Arabidopsis were still much lower than these in unstressed Haberlea. While 
all species accumulated sucrose during chilling, that accumulation was transient in H. 
rhodopensis and A. thaliana but sustained in T. halophila after the return to optimal 
temperature. Thus, Haberlea's metabolome appeared primed for chilling stress but 
the low temperature acclimation induced additional stress-protective mechanisms. A 
diverse array of sugars, organic acids, and polyols constitute Haberlea's main metabolic 
defence mechanisms against chilling, while accumulation of amino acids and amino acid 
derivatives contribute to the low temperature acclimation in Arabidopsis and Thellungiella. 
Collectively, these results show inherent differences in the metabolomes under the 
ambient temperature and the strategies to respond to low temperature in the three 
species. 

Keywords: Arabidopsis thaliana, Haberlea rhodopensis, low temperature stress, metabolite profiling, Thellungiella 
halophila 



INTRODUCTION 

The small but diverse group of resurrection plants exhibit a 
remarkable adaptation to extreme drought stress. In absence of 
water supply, they can tolerate desiccation of their vegetative 
tissues to air dried state and quickly regain normal appear- 
ance and metabolism upon rehydration (Dinakar et al., 2012; 
Gechev et al., 2012). Haberlea rhodopensis is a desiccation-tolerant 
species, perennial herbaceous plant endemic to several mountains 
in the Balkan Peninsula in South-Eastern Europe (Gechev et al., 
2013a). It is also an ancient plant, a glacial relic, which might 
have acquired its defence mechanisms a long time ago. As it is 
exposed to the harsh winter conditions and subzero temperatures 



in these latitudes, this species additionally evolved mechanisms to 
withstand chilling and freezing stress. 

Earlier studies on resurrection plants indicated that com- 
plex and diverse mechanisms can contribute to their desicca- 
tion tolerance. These include alterations of sugar metabolism, 
reconfiguration of the cell wall, inhibition of growth and pho- 
tosynthesis, rapid induction of late embryogenesis abundant 
(LEA) and small heat shock proteins, accumulation of pheno- 
lic antioxidants, upregulation of antioxidant enzymes, aldehyde 
dehydrogenases, and other protective enzymes (Kirch et al., 2001; 
Mowla et al, 2002; Battaglia et al., 2008; Rodriguez et al, 2010; 
Van Den Dries et al., 2011; Moore et al, 2012; Gechev et al., 
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2013a). Transcriptional re-programming and metabolome re- 
adjustments are important elements of this stress defence strategy 
(Rodriguez et al, 2010; Oliver et al, 2011; Yobi et al., 2012, 2013; 
Gechev et al., 2013a). However, little is known about the molecu- 
lar responses of resurrection species to low temperatures and no 
resurrection species has been investigated in terms of metabolome 
reconfiguration during low temperature stress. 

Exposure to freezing environments leads to serious damage 
of the plant cell by ice formation and dysfunction of cellular 
membranes. Many plant species increase freezing tolerance dur- 
ing exposure to non-freezing low temperature by a process known 
as "cold acclimation." The molecular basis of this process has 
been extensively studied in Ambidopsis thaliana, which is con- 
sidered sensitive to cold stress, and the contribution of particular 
metabolites including compatible solutes and the transcriptional 
regulatory network has been elucidated. For instance, accumu- 
lation of sugars is considered to play an important role in cold 
acclimation (Hannah et al., 2006) and transcription factors such 
as CBF3/DREB1A play a central role to control this process 
(Cook et al., 2004). Thellungiella halophyla is a close relative of 
A. thaliana that has been suggested to possess the characteris- 
tics of an extremophile, i.e., high tolerance to salinity, freezing, 
nitrogen-deficiency, and drought stress (Lee et al., 2012). For 
this reason, Thellungiella has been analyzed in comparison to 
Arabidopsis to elucidate the mechanisms that confer tolerance 
against abiotic stress. Although some accessions of Thellungiella 
are not extremophile with regard to freezing tolerance, others, 
including Yukon, show significantly higher tolerance than any 
Arabidopsis accessions (Lee et al, 2012). The metabolite profiling 
data show different metabolic adaptation strategies between these 
two species (Lee et al., 2012), indicating specific cold acclima- 
tion processes which lead to the different levels of cold tolerance. 
Recent studies on Picea sitchensis and Fragaria vesca confirmed 
the notion that specific cold acclimation processes exist (Dauwe 
et al., 2012; Rohloff et al., 2012). The desiccation tolerance of 
H. rhodopensis outperforms both Thellungiella and Arabidopsis. 
Furthermore, Haberlea can withstand freezing temperatures, sug- 
gesting distinctive cold acclimation strategies allowing high freez- 
ing tolerance in this species. The main aim of this study was 
to reveal the metabolic changes of H. rhodopensis during low 
temperature treatment and subsequent return to optimal growth 
temperature. Comparison of the strategies for metabolic adap- 
tation to cold in H. rhodopensis, T. halophila and A. thaliana as 
representatives of resurrection plants, extremophiles and non- 
extremophiles, respectively, was carried out to highlight the dif- 
ferences and the common pathways these species employ to 
adapt to low temperatures. The results suggest the importance 
of metabolite composition under non-stress conditions as a pre- 
adaptation strategy and point out the diverse low-temperature 
stress responses in these three species which likely contribute to 
the different levels of stress tolerance. 

MATERIALS AND METHODS 

PLANT MATERIAL, GROWTH CONDITIONS, AND LOW TEMPERATURE 
TREATMENT 

A. thaliana ecotype Col-0 was obtained from the Nottingham 
Arabidopsis Stock Centre (NASC, http://arabidopsis.info/); H. 
rhodopensis was initially collected from the Rhodope mountains 



and subsequently maintained in a climate-controlled room on 
soil taken from its natural habitat as described (Gechev et al., 
2013a,b); T. halophila ecotype Yukon was obtained from Dr. 
Yang-Ping Lee and Dr. Dirk Hintcha, Max-Planck Institute of 
Molecular Plant Physiology, Potsdam-Golm, Germany. 

Plants were grown in a climate room on soil at 21°C, 40 u,mol 
mT 2 s _1 light intensity, 16/8 light/dark photoperiod, and rel- 
ative humidity 70%. Rosette leaves from all three species were 
used as samples. Low temperature stress was applied by placing 
the plants in a plant growth chamber at 4°C, 40 u,mol m~ 2 s 
light intensity, 16/8 light/dark photoperiod, and relative humid- 
ity 70%. Samples were taken after 3 days of chilling. Subsequently, 
plants were transferred back to the normal growth conditions 
(21°C, 40 ixmol. m~ 2 s _1 light intensity, 16/8 light/dark photope- 
riod, and relative humidity 70%), and samples were taken after 3 
days of recovery. The duration of the low temperature treatment 
and the recovery period were chosen because longer stress periods 
interfered with the development of Arabidopsis and are supposed 
to induce secondary effects on the metabolite profile. In all cases, 
plant material was immediately frozen in liquid nitrogen, ground 
into fine powder and stored at — 80° C until further analysis. Six 
biological replicates were used for the analyses. 

MEASUREMENTS OF ANTHOCYANINS, MALONDIALDEHYDE, AND 
REDUCED AND OXIDIZED GLUTATHIONE 

The anthocyanins and malondialdehyde were measured 
photometrically as described by Gechev et al. (2013a,b). 
Briefly, anthocyanins were extracted with 1% HC1 in methanol 
and the anthocyanin content was determined by reading the 
absorbance at 530 nm. Correction for non-specific absorption 
of photosynthetic pigments was done at 657 nm and the final 
anthocyanin amount, calculated as A530 — 0.25 A^s7, was nor- 
malized by the fresh weight of the samples. Malondialdehyde 
was extracted using 1 ml 0.25% thiobarbituric acid dissolved 
in 10% trichloroacetic acid. After heating the extracts at 85°C 
for 30min and rapid chilling on ice, the pellets were removed 
by centrifugation and the specific absorbance read at 532 nm 
(the peak of malondialdehyde-thiobarbituric acid complex). 
Correction for non-specific absorbance at 600 nm was made, 
malondialdehyde concentration calculated using an extinction 
coefficient E532 155 mM -1 cm -1 , and values normalized by the 
fresh weight of the samples. 

Glutathione, total and oxidized, was measured by an enzy- 
matic assay essentially as described by Mehterov et al. (2012). 
The method relies on the GR-dependent reduction of DTNB, 
monitored at 412 nm. Briefly, samples were homogenized in 
1 ml 5% sulfosalicylic acid (Sigma-Aldrich, St. Louis, Missouri, 
USA) made in 0.1 M potassium phosphate buffer (pH 7.6/5 mM 
EDTA) on ice. Aliquots of neutralized extract were mixed with 
1.2 mM DNTB (Sigma-Aldrich, St. Louis, Missouri, USA) and 
0.3 mM NADPH in 0.1 M potassium phosphate buffer (pH 
7.6/5 mM EDTA) and the reaction was started by the addition 
of 1 U glutathione reductase (Sigma-Aldrich, St. Louis, Missouri, 
USA). The increase in A412 was monitored for 1 min. Oxidized 
glutathione was measured by the same principle after incubation 
of neutralized extract with 2 |xl 2-vinylpyridine (Sigma-Aldrich, 
St. Louis, Missouri, USA) for 1 h at room temperature to com- 
plex the reduced glutathione. To remove excess 2-vinylpyridine, 
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the derivatized solution was treated with diethyl ether. Reduced 
glutathione was determined as the difference between total and 
oxidized glutathione. 

RNA ISOLATION AND QUANTITATIVE RT-PCR 

RNA from frozen leaf material was extracted with Trizol Reagent 
(Invitrogen, Life Technologies, Carlsbad, California, USA) 
according to the manufacturer's recommendations. Ten micro- 
grams of total RNA was treated with DNA-free™ Kit (Ambion) 
to remove any DNA contamination. RNA integrity was checked 
on 1% (w/v) agarose gel and concentration measured with a 
Nanodrop ND-2000 Spectrophotometer before and after DNAse 
I digestion. Additionally, the quality and integrity of the RNA 
samples were analyzed on an RNA 6000 Lab-on-a-Chip using the 
Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). 
cDNA was synthesized from 2 u,g of total RNA using RevertAid™ 
First Strand cDNA Synthesis Kit (Fermentas, Thermo Fisher 
Scientific, Waltham, Massachusetts, USA) with oligo-dT primers, 
according to the manufacturer's instructions. 

For qRT-PCR analysis, A. thaliana and T. halophila chilling- 
and drought stress-responsive genes were selected based on lit- 
erature data and information from colleagues (Carvallo et al., 
2011; Zou et al, 2011; Y. Ping Lee, personal communica- 
tion). These included genes encoding two drought stress tran- 
scription factors DREB2A and DREB2B in Arabidopsis and 
DREB2B in Thellungiella, genes encoding LEA proteins in 
both species, the cold stress-inducible COR15A, COR47, RD29A 
genes in Arabidopsis, and ELIP2, ADH1, and RD29B genes 
in Thellungiella. Some of these genes, such as RD29A/RD29B 
(responsive to dehydration), are reported to respond to dehydra- 
tion as well as to low temperature stress, due to the cellular dehy- 
dration that occurs at low temperatures, while others are believed 
to be more stress-specific. As there were no studies on the molec- 
ular responses to low temperatures in Haberlea, genes for the 
transcriptional analysis were chosen based on our previous anal- 
ysis during dehydration. The genes were selected from those most 
regulated by dehydration (Gechev et al, 2013a) and homologs of 
cold-inducible genes in other species including the temperature- 
induced lipocalin and the cold-induced glucosyl transferase. The 
genes of Arabidopsis and Thellungiella are very close to each other 
in terms of sequence homology, as expected for closely related 
species. The Habeiiea's genes are very divergent from the other 
two species but nevertheless serve as useful stress markers. 

Quantitative real-time PCR (qRT-PCR) analysis was per- 
formed using an ABI PRISM 7900 HT PCR instrument 
(Applied Biosystems, Life Technologies, Carlsbad, California, 
USA). Primers for the qRT-PCR analysis were designed using 
the Primer3 software. The genes and corresponding primer 
pairs from the three species used for the analysis are listed in 
Supplemental Table 1. All reactions contained 10 (xL of SYBR 
Green Master Mix (Applied Biosystems, Life Technologies, 
Carlsbad, California, USA), 25 ng of cDNA, and 200 nM of each 
gene-specific primer in a final volume of 20 (xL. The qRT-PCR 
reactions were executed using the following program: 50°C for 
2min, 95° C for lOmin, followed by 40 cycles of 95°C for 15 s 
and 60° C for 1 min. 

Relative mRNA abundance was calculated using the com- 
parative 2~ AACt method and normalized to the corresponding 



reference gene levels (Schmittgen and Livak, 2008). Fold changes 
in gene expression were calculated and the resulting data sets 
were log transformed and visualized by Multi Experiment Viewer 
(MEV)-created heat maps. Two biological and two technical 
repetitions were performed for each gene. 

METABOLITE PROFILING 

Profiling of primary metabolites was conducted using an 
established gas chromatography mass spectrometry (GC-MS) 
protocol exactly as described by Lisec et al. (2006). The 
metabolites were extracted from 100 mg frozen material by 
methanol and polar metabolites were isolated by phase sep- 
aration using chloroform. The polar phase was taken and 
aliquots of 150 |xl were dried for further analysis. The metabolites 
were derivatized with methoxyamine-HCl and N-Methyl-N- 
(trimethylsilyl)trifluoroacetamide and subjected to GC-MS anal- 
ysis. Samples were analyzed by a Gas chromatograph, 6890N 
(Agilent Technologies, Santa Clara, CA) connected to Pegasus III 
time -of- flight mass spectrometer (Leco Instruments, St. Joseph, 
MI) using a MDN-35 capillary column (Macherey-Nagel, Diiren, 
Germany). Chromatograms and mass spectra were evaluated 
by Chroma TOF®4.2 (Leco, St Joseph, MI) and TagFinder 4.0 
(Luedemann et al, 2008) for the quantification and annota- 
tion of the peaks using the MPI Golm Metabolome Database 
(GMD, http://gmd.mpimp-golm.mpg.de/, Kopka et al., 2005). 
The parameters used for the identification of the metabolite are 
summarized following the way recommended in Fernie et al. 
(2011) as Supplemental Table 3. The amount of metabolites was 
analyzed as relative metabolite abundance calculated by normal- 
ization of signal intensity to that of ribitol which was added as 
an internal standard and then by the fresh weight of the mate- 
rial. The whole dataset is provided in Supplemental Table 4. 
Metaboanalyst (www.metaboanalyst.ca, Xia et al., 2012) was used 
for data analysis including principal component analysis (PCA). 

RESULTS 

CHARACTERIZATION OF THE PHYSIOLOGICAL RESPONSES OF 
H. rhodopensis, T. halophyla, AND A. thaliana TO CHILLING 
TREATMENT AND SUBSEQUENT RECOVERY 

To evaluate the influence of low temperatures on the three 
species, plants were inspected for any visible damage and a 
number of physiological parameters were measured: malondi- 
aldehyde levels, which are indicators of lipid peroxidation and 
oxidative stress; chlorophyll pigments, which normally decrease 
during severe stress; reduced and oxidized glutathione, which 
increase as a result of various stresses and their ratio indicates 
the redox status of the cell. Exposure of H. rhodopensis, T. halo- 
phyla, and A. thaliana to 4°C for 3 days did not cause any visible 
tissue damage or cell death (data not shown), nor any signifi- 
cant increase in lipid peroxidation as judged by determination 
of malondialdehyde levels (Figure 1). Furthermore, no visible 
decrease in turgor or wilting was observed. In control condi- 
tions, Haberlea and Thellungiella had 2-fold higher levels of glu- 
tathione than Arabidopsis (Figure 1). While the reduced (GSH) 
and oxidized (GSSG) glutathione remained unchanged dur- 
ing cold treatment and recovery in Haberlea and Thellungiella, 
both GSH and GSSH increased significantly in Arabidopsis 
during cold. In Arabidopsis, GSH then returned to initial 
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FIGURE 1 | Malondialdehyde (MDA), anthocyanins, reduced 
glutathione (GSH), and oxidized glutathione (GSSG) levels in 
A. thaliana, H. rhodopensis, and T. halophila exposed to low 
temperature stress and subsequent recovery. Plants grown at optimal 
conditions (controls) at 21°C, 40iimol. m~ 2 s~ 2 light intensity, 16/8 
light/dark photoperiod, and relative humidity 70%, were subjected to 
chilling (4°C) for 3 days and then returned to 21 °C for 3 days for 



recovery. (A), MDA content; (B), anthocyanins; (C), GSH levels; (D), 
GSSG levels. MDA, GSH, and GSSG are expressed as nmol per gram 
fresh weight (FW), while anthocyanins are presented as the difference 
in absorbencies between 530 and 567 nm (A 5 3o — 0.25 A 5 67) per 100 mg 
fresh weight. Asterisks indicate significant difference between controls 
and low temperature-treated or recovered plants (P < 0.05, student's 
f-test). The data are means ± s.e.m. of three replicates. 



values on recovery, while GSSH decreased but was still higher 
than in unstressed controls (Figure 1). Unstressed Haberlea and 
Thellungiella also displayed 3- and 4-fold higher levels of antho- 
cyanins than Arabidopsis, respectively (Figure 1). In Arabidopsis 
and Thellungiella, anthocyanins increased on recovery (Figure 1). 

LOW TEMPERATURE- AND DEHYDRATION-RESPONSIVE MARKER 
GENES ARE INDUCED BY CHILLING TREATMENT IN H. rhodopensis, A. 
thaliana, AND 7! halophyla 

Cold treatment of H. rhodopensis, T. halophyla, and A. thaliana 
resulted in induction of low temperature- and dehydration- 
responsive marker genes in all three species (Figure 2, 
Supplemental Table 2). The data shows that cold treatment 
at 4°C for 3 days induced significantly the COR15A, COR47, 
and RD29A genes of Arabidopsis and RD29B of Thellungiella, 
as well as the LEA gene of Arabidopsis and the ELIP2 and ADH1 
genes of Thellungiella (Figure 2, Supplemental Table 2). The 
DREB2A and DREB2B genes of Arabidopsis were not upreg- 
ulated, indicating that the low temperature treatment did not 
cause dehydration. In Haberlea, four of the six selected genes 
were induced by the low temperature, including the temperature- 
induced lipocalin and the cold-induced glucosyl transferase. A 
protein phosphatase gene had the highest level of induction; 
interestingly, the same gene showed the strongest induction by 
dehydration (Gechev et al., 2013a). With the exception of this 
protein phosphatase, the levels of all Arabidopsis, Thellungiella, 
and Haberlea genes returned to normal values upon recovery 



from chilling (Figure 2). Taken together, these results indicate 
all three species induced cold response by the low-temperature 
treatment and recovered following the treatment. 

METABOLITE PROFILES OF H. rhodopensis, A. thaliana, AND T. 
halophyla UNDER NON-STRESS CONDITIONS AND THEIR 
RECONFIGURATIONS DURING LOW TEMPERATURE AND SUBSEQUENT 
RECOVERY 

The measurement of the relative metabolite levels in H. 
rhodopensis and subsequent comparison with A. thaliana and 
T. halophyla showed substantial differences between the three 
species. PCA revealed the differences in the metabolite pro- 
files of the three species in normal growth condition rather 
than the changes during cold treatment in the same species 
(Figure 3). A. thaliana displayed global changes in metabolite 
levels during the cold treatment while the other two species did 
not (Figure 3). The metabolite levels under control conditions 
were compared in Figure 4. A. thaliana had much higher levels 
of putrescine and fumarate than the other two species. Unique 
features of the Haberlea metabolome were the high level of 
many sugars and sugar alcohols including glucose, fructose, 
sucrose, trehalose, rhamnose, raffinose, galactinol, myo-inositol, 
and sorbitol (Figure 4). This plant also accumulated some 
organic acids including two TCA cycle intermediates, citrate, and 
2-oxoglutarate. Thellungiella, on the other hand, had the highest 
levels of amino acids including arginine, asparagine, threonine, 
pyroglutamate, histidine, phenylalanine, valine, glutamine, 
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FIGURE 2 | Drought- and cold-responsive genes are induced by the low 
temperature treatment in A. thaliana, H. rhodopensis, and T. halophila. 



Plants grown at optimal conditions at 21 °C, 40 ixmol. m 



light 



intensity, 16/8 light/dark photoperiod, and relative humidity 70%, were 
subjected to chilling (4°C) for 3 days and then returned to 21 °C for recovery. 
Fully hydrated plants grown at 21 °C were used as a reference point for the 
qRT-PCT analysis. Red colors and positive values depict induced genes 
while green colors and negative values depict repressed genes. The data 
are log ratios, average means of two replicates. The genes and primer 
combinations used are given in Supplemental Table 1, while the exact 
values of induction/repression are given in Supplemental Table 2. 



lysine, isoleucine, ornithine, tyrosine, beta-alanine, serine, 
tryptophane, and proline (Figure 4). Another observation was 
the lower levels of ascorbate and dehydroascorbate in Haberlea, 
than in Arabidopsis and Thellungiella (Figure 4). 

In addition to these intrinsic metabolic differences, the three 
species showed distinct metabolic responses to the low tem- 
perature treatment. All of them accumulated sucrose, fructose 
and glucose during cold treatment indicating that accumu- 
lation of these sugars is a common response in cold accli- 
mation. Raffinose was detected only after cold treatment in 
Arabidopsis and Thellungiella, while it was abundant already 
under normal conditions in Haberlea (Supplemental Table 4). 
Whilst maltose was not detectable in Arabidopsis, it was clearly 
detectable in Thellungiella and Haberlea samples following 
cold treatment and recovery, where it accumulated to much 
higher levels (Supplemental Table 4). In addition, aspartate tran- 
siently accumulated under cold treatment in all three species 
although the relevance of this to cold acclimation is currently 
unclear. As suggested by PCA analysis (Figure 3), Arabidopsis 
transiently accumulated many metabolites including galactinol 
and proline (Figure 5). An increase in alanine, putrescine and 
pyruvate was observed both in Arabidopsis and Thellungiella 
(Figure 5). In many cases, the metabolite changes measured 
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FIGURE 3 | Principal component analysis of metabolite profiles of 
leaves from A. thaliana (A), T. halophyla (T) and H. rhodopensis (H) 
under optimal growth condition (0) and following cold treatment (1) 
and subsequent recovery (2). Each point represents an individual 
biological replicate. Plotting of the first and second component is shown. 
The circles indicate the 95% confident regions. 



during low temperature treatment were transient. In contrast, 
there were several examples of sustained accumulation includ- 
ing sucrose, proline, urea and 4-hydroxyproline in Thellungiella 
(Figure 5). 

DISCUSSION 

Despite the numerous studies on resurrection species in the past 
decade, there is no detailed information on the response of res- 
urrection species to other types of abiotic stresses except drought. 
In particular, no information is available on the molecular mech- 
anisms of low temperature tolerance in any of the angiosperm 
resurrection species. We chose to investigate H. rhodopensis, as 
this perennial plant can withstand harsh winters in its native 
habitat and a number of genes associated with low tempera- 
ture responses are induced by dehydration, suggesting a possible 
cross-protection (Gechev et al., 2013a). While our main goal was 
to investigate the metabolic reconfiguration of Haberlea during 
low temperature treatment and subsequent recovery, compara- 
tive analysis with Arabidopsis and Thellungiella highlighted the 
species-specific metabolic responses as well as the protective 
mechanisms conserved among evolutionary distant organisms. 
Due to the differences in morphological and physiological prop- 
erties of these species, direct comparison of metabolite levels 
should be done with special care. However, prominent differences 
(considerered as qualitative) in the metabolic levels and the accu- 
mulation pattern during the time course should reflect metabolic 
feature of the species. 
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FIGURE 4 | Comparison of metabolite levels in A. thaliana (blue), 77 halophyla (red) and H. rhodopensis (green) under pre-stress condition. Metabolite 
levels are normalized by the means of all samples and presented as mean + s.e.m. of six biological replicates. 
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FIGURE 5 | Changes of metabolite levels during cold treatment and 
subsequent recovery. The levels of metabolites were normalized by those at 
control conditions in each plant species. Ten metabolites which showed 
statistically significant changes following the treatment are shown. The 



values are presented as means ± s.e.m. of six biological replicates. Asterisks 
indicate significant differences by student's f-test (P < 0.05) of the treated 
samples compared to the control time points in each species. Blue diamond, 
A. thaliana; red square, H. rhodopensis; green triangle, T halophyla. 



LOW TEMPERATURE CAUSES COLD ACCLIMATION RESPONSES BUT 
NO SEVERE OXIDATIVE STRESS IN H. rhodopensis, A. thaliana, AND T. 
halophyla 

The duration and magnitude of the low temperature treatment 
caused neither severe oxidative stress, as judged by the unal- 
tered malondialdehyde and chlorophyll levels, nor cell death in 



any of the three species. Yet, the low temperatures imposed 
prominent changes in gene expression and metabolite levels of 
Haberlea, Arabidopsis, and Thellungiella, indicating that acclima- 
tion responses took place in all three species. For Haberlea, these 
results are original. Although a number of studies investigated dif- 
ferent aspects of low temperature responses of Arabidopsis and 
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Thellungiella, no single study covered all aspects of low tempera- 
ture stress (Cook et al., 2004; Hannah et al., 2006; Alcazar et al., 
2011; Carvallo et al., 2011; Lee et al, 2012; Zuther et al, 2012). 
Furthermore, even small differences between growth conditions 
or/and experimental set-ups may result in notable differences in 
metabolite levels, hence the two species were also included in 
the experiments, instead of retrieving data for them from the 
literature. 

The higher basal levels of anthocyanins in Haberlea and 
Thellungiella compared with Arabidopsis (3- and 4-fold higher, 
respectively) may be a part of a pre-adaptation strategy, ensur- 
ing a higher level of basal protection, although interpretation 
of such cross-species comparisons need to be done carefully. 
Anthocyanins are well-known antioxidants and protectors against 
oxidative stress. Recently, tomato plants expressing higher levels 
of anthocyanins have been demonstrated to produce fruits with 
expended shelf life, delayed ripening, and increased resistance to 
the fungal pathogen Botrytis cinerea, due to the altered spreading 
of ROS burst after infection (Zhang et al, 2013). Furthermore, 
the oxidative stress-tolerant Arabidopsis mutant oxrl was shown 
to have much higher levels of anthocyanins than the wild type and 
much stronger induction of anthocyanin synthesis upon oxidative 
stress (Gechev et al, 2013b). 

The higher levels of the prominent antioxidant and redox 
regulator glutathione in Haberlea and Thellungiella may also con- 
tribute to the higher basal stress tolerance of the two species 
relative to Arabidopsis. It is interesting to note that the lev- 
els of reduced glutathione remain unchanged in Haberlea and 
Thellungiella during chilling and subsequent recovery while the 
reduced glutathione increases 2-fold in low temperature treated 
Arabidopsis to reach the levels of the other two species, but then 
drops again during recovery. Furthermore, oxidized glutathione 
levels do not change in Haberlea during chilling and recovery but 
increased 3-fold when Arabidopsis was subjected to chilling. This 
collectively implies that the low temperature changes the oxida- 
tive status of the more sensitive Arabidopsis and has less effect on 
the more tolerant Thellungiella and Haberlea. 

LOW TEMPERATURE- AND DEHYDRATION-RESPONSIVE MARKER 
GENES ARE INDUCED BY CHILLING IN H. rhodopensis, A. thaliana, AND 
7! halophyla 

Despite Haberlea's potential in plant biotechnology and biomed- 
ical science, the genome of this species remains unsequenced 
and very limited information is currently available concerning its 
genes (Apostolova et al., 2012; Georgieva et al., 2012). Recently, 
the first comprehensive transcriptome and metabolite profiling 
study of H. rhodopensis indicated the types of genes present in 
Haberlea and how they are expressed during dehydration and 
rehydration (Gechev et al, 2013a). The low temperature- and 
drought-responses of Arabidopsis and Thellungiella are well stud- 
ied (Cook et al, 2004; Hannah et al, 2006; Thomashow, 2010; 
Carvallo et al., 2011; Zou et al, 2011; Lee et al, 2012), which 
allowed us to select proper marker genes for use as indicators 
of stress. The results from this study provide valuable infor- 
mation concerning low temperature-inducible genes that can 
be used as markers for abiotic stress in Haberlea. The cold- 
inducible glucosyl transferase (Figure 2, Supplemental Table 2) 
seems to be specific for low temperature as our expression analysis 



shows that it is induced under low temperature but not signifi- 
cantly by drought (Gechev et al., 2013a). The other three genes 
that were upregulated after chilling, those encoding a protein 
phosphatase, a MADS box protein, and a temperature-induced 
lipocalin, are also responsive to drought and desiccation and can 
be used as general markers for these two types of abiotic stress. 
Further expression studies on these genes during other types 
of stress will show if they can be utilized as general markers 
of abiotic stress or they respond only to low temperature- and 
drought-induced osmotic stress. The protein phosphatase gene is 
of particular interest, as it was the most strongly induced gene by 
both low temperature and drought stress (Figure 2; Gechev et al, 
2013a). 

METABOLIC CHANGES INDICATED DISTINCTIVE STRATEGIES OF COLD 
ACCLIMATION IN H. rhodopensis, A. thaliana, AND 7! halophyla 

The three species have very distinct metabolic profiles, clearly vis- 
ible already in the absence of stress. This suggests the importance 
of basal metabolic composition for cold tolerance. Metabolic pre- 
adaptation is considered as a major factor affecting the tolerance 
to stresses in species and cultivars (Sanchez et al., 2011). 

Sugar metabolism in particular plays a major role in several 
types of abiotic stress responses, especially in tolerance against 
drought, osmotic stress, and chilling. Haberlea was revealed to 
possess unique sugar metabolism with a variety of sugars and 
sugar alcohols of much higher levels than the other two species. 
Thellungiella, on the other hand, accumulates more sugars and 
sugar alcohols than Arabidopsis (Gong et al., 2005). Higher 
expression levels of stress-tolerant genes in Thellungiella and 
accumulation of several compounds that have protective func- 
tions in the presence of osmotic imbalance are evident even 
under pre-stress conditions (Gong et al, 2005). Recent works in 
Arabidopsis suggested that it may not be a specific sugar that is 
important for plant freezing tolerance, but rather that sugars in 
general may constitute a highly redundant cryoprotective system 
(Korn et al., 2008; Zuther et al., 2012). The accumulation of a 
large variety of metabolites is likely to contribute to establish- 
ment of a robust system to cope with environmental stresses in 
Haberlea. Particularly, the higher levels of maltose indicate more 
intensive starch breakdown. This, together with the much higher 
levels of the monosaccharides glucose and fructose, appears to 
fuel sucrose synthesis in Haberlea. In the other two species, the 
amount of fructose and glucose is much lower, although increases 
are clearly visible during low temperature conditions. Sucrose, 
fructose and glucose accumulation seems to be a general response 
of all the three species to low temperature stress and to osmotic 
stress in particular (Gechev et al., 2013a). Sucrose accumulation is 
known as a common response to various environmental stresses 
(Obata and Fernie, 2012). The much higher levels of sucrose and 
trehalose in Haberlea than in the other two species imply a role 
of these two non-reducing disaccharides in both drought and 
low temperature stress tolerance. They are known to accumu- 
late in resurrection plants during dehydration (Drennan et al., 
1993; Ingram and Bartels, 1996; Norwood et al, 2000, 2003; 
Martinelli, 2008) and can serve as osmoprotectants of biological 
membranes and can stabilize macromolecular structures (Crowe 
et al., 1992; Martinelli, 2008). Besides this, both sugars have also 
signaling properties in lower concentrations. Trehalose and its 
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derivative, trehalose-6-phosphate, are central metabolic regula- 
tors in Arabidopsis, influencing carbohydrate status, growth, and 
energy metabolism (Schluepmann et al., 2003; Lunn et al., 2006; 
Smeekens et al., 2010). Trehalose-6-phosphate stimulates ADP- 
glucose pyrophosphorylase, promoting starch synthesis, while 
trehalose has the opposite effect stimulating starch breakdown 
(Schluepmann et al., 2003; Smeekens et al., 2010). Since the 
sugars should be highly accumulated to function as a compati- 
ble solute, the qualitatively different levels of trehalose suggests 
its different roles in the three species, namely as a compatible 
solute in Haberlea and as a signaling molecule in Arabidopsis and 
Thellungiella as proposed elsewhere (Carillo et al., 2013). 

Haberlea has also much higher levels of raffmose, myo- 
inositol and galactinol, which are precursors of raffmose fam- 
ily oligosaccharides. The higher amount of galactinol observed 
in Haberlea may be directly linked to protection against abi- 
otic stress. Galactinol and raffmose can protect against oxidative 
stress, which is a consequence of many abiotic stresses, includ- 
ing chilling (Nishizawa et al., 2008). In Arabidopsis, oxidative 
stress-mediated induction of galactinol synthases and raffmose 
synthases is governed by the heat shock transcription factor 
HsfA2 and eventually leads to elevated levels of galactinol and 
raffmose (Nishizawa et al., 2008). In a recent study, high lev- 
els of galactinol as well as expression of dehydrin and alcohol 
dehydrogenase genes during cold acclimation correlated with 
low temperature tolerance in F. vesca (Davik et al., 2013). The 
raffmose family oligosaccharides, including raffmose, stachyose, 
and verbascose, have been shown to accumulate during osmotic 
stresses such as drought and chilling and protect from these 
stresses by water replacement and vitrification (Norwood et al., 
2003; Farrant et al., 2007). Indeed, the high levels of raffmose 
accumulated in Arabidopsis during chilling may be the manner 
in which this species defends itself from low temperature stress. 
Haberlea's pathway may also be shifted toward synthesis of the 
more complex stachyose and verbascose, as it is during drought 
and desiccation (Gechev et al., 2013a). 

Induction of putrescine by the low temperature treatment, 
observed in Arabidopsis and Thellungiella but not in Haberlea, 
may be a part of the metabolic reconfigurations of the two 
species to counteract the stress. Recently, induction of putrescine 
by cold acclimation was observed in F. vesca (Rohloff et al., 
2012). They also found cold-induced accumulation of aspar- 
tate, which is observed in our case in all three species. Levels 
of polyamines including putrescine raise under cold exposure in 
many species. There is considerable evidence for the important 
roles of polyamine including putrescine in plant defence against 
cold and other abiotic stresses although the mode of action is still 
elusive (Alcazar et al., 2011). 

Unique features of Thellungiella metabolome were the very 
high levels of amino acids. Some amino acids are known to con- 
tribute to the tolerance to abiotic stresses. Proline is one of the 
well-documented osmoprotectants and some other amino acids 
including branched chain amino acids are accumulated under 
various stress conditions (Obata and Fernie, 2012). Interestingly, 
Thellungiella further accumulated proline and hydroxyproline 
during recovery from cold, which may be relevant for protection 
against subsequent or sustained low temperature stress. Proline 
is a well-known osmoprotectant that accumulates in response to 



drought or chilling in a number of species but Haberlea does 
not utilize this amino acid for protective purposes, as proline 
levels in Haberlea are low and, unlike Arabidopsis, do not ele- 
vate during chilling. These observations suggest that Arabidopsis 
and Thellungiella take advantage of amino acids and amino acid 
derivertives for stress adaptation in contrast to Haberlea, which 
accumulates sugars. The reason of this preference is unclear since 
the known functions of sugars and amino acids in stress tol- 
erance are similar (Obata and Fernie, 2012). It may be due to 
the energetic and/or nutritional cost to produce large amount of 
metabolites belonging to these groups. 

COMPARATIVE ANALYSIS BETWEEN LOW TEMPERATURE- AND 
DEHYDRATION-INDUCED METABOLIC RESPONSES OF Haberlea 
rhodopensis 

The changes in metabolite levels of Haberlea during low tem- 
perature stress were compared with the recently conducted 
metabolome analysis of this species during drought stress and 
desiccation (Gechev et al., 2013a). Low temperature and drought 
are distinct abiotic stresses but both of them lead to loss of cellu- 
lar water, hence specific as well as common metabolic responses 
are expected. The most obvious common response during both 
stresses was the massive accumulation of sucrose, which identifies 
this metabolite as a common protector against the two different 
types of osmotic shock in Haberlea. In both cases, sucrose accu- 
mulation is accompanied by elevated maltose levels, implying that 
starch degradation occurs during drought and low temperature 
stress — most probably as a carbon and energy source. During 
drought, this is accompanied by reduction of glucose and sucrose 
levels. During cold and especially upon recovery, however, the two 
sugars actually increase, especially upon recovery. 

Apart from these few examples, most metabolites responded to 
the two stresses differently, suggesting very specific metabolome 
reconfigurations for drought and low temperature overall. Many 
key metabolites like trehalose, proline, citrate, succinate, and 
malate remained unchanged during low temperature, while all of 
them decreased during drought and desiccation. These metabo- 
lites, some of them well-known osmoprotectors and/or signaling 
molecules, may not be involved in acquiring drought tolerance 
but could well be part of the cold stress defence. Asparagine 
and aspartate levels decreased during drought/desiccation but 
increased during low temperature treatment. Aspartic acid was 
also induced by cold in F. vesca (Rohloff et al., 2012) and could 
contribute to the low temperature defence in H. rhodopensis. 
Proline and several organic acids related to the tricarboxylic acid 
cycle, including citrate, succinate, and malate, decreased during 
drought and desiccation but remained constant during low tem- 
perature stress. Trehalose levels also remained constant during 
chilling, while they dropped during dehydration and subsequent 
rehydration. GABA, a well-known stress signaling molecule and 
growth regulator, may on the other hand be involved primarily 
in the molecular mechanisms of desiccation tolerance, as its lev- 
els dramatically increase during desiccation but not during low 
temperature stress. 

In conclusion, the three species have intrinsic differences in 
their metabolomes in the absence of stress and respond differently 
to chilling, implying unique strategies to counteract low tem- 
perature stress. Haberlea's high levels of galactinol, myo-inositol, 
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sorbitol, and many sugars provide this species with steady-state 
metabolome which is configured to encounter the consequences 
of the low temperature stress already during the optimal growth 
conditions. While the oxidative stress-protective properties of 
galactinol and the osmoprotective as well as signaling properties 
of sucrose and trehalose are well-documented, the function of the 
other sugars such as fucose, rhamnose, and other monosaccha- 
rides in low temperature stress remains to be studied. Glucose 
and fructose, for example, may have dual functions: on one side, 
they are substrate for the increased sucrose synthesis, but on 
the other side they together with the rhamnose and fucose may 
play a role in the reconfiguration of the cell wall polysaccha- 
rides. Arabidopsis could utilize another defensive strategy based 
on transient cold-induced accumulation of sucrose, putrescine, 
and proline. It should also be noted that sugars in general may be 
involved in the cold acclimation rather than a specific sugar (Korn 
et al., 2008; Zuther et al, 2012). The high number of uniden- 
tified metabolites exclusively present in H. rhodopensis (Gechev 
et al., 2013a) suggests that there may be unique compounds 
this species synthesizes to protect itself from abiotic and oxida- 
tive stress. The accumulation of a large variety of metabolites is 
likely to contribute to the establishment of a robust system to 
cope with environmental stresses in Haberlea. Further study of 
the chemical identity of the unidentified metabolites may thus 
reveal new compounds with powerful stress protective functions 
(Gechev et al., 2013a). Thellungiella utilizes a third strategy, based 
on pre-adaptation using amino acids and polymines, transient 
cold-induced accumulation of putrescine and amino acids such as 
alanine and aspartate, and sustained accumulation of sucrose and 
proline. The latter may give Thellungiella an advantage relative to 
Arabidopsis in future stress encounters. 
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